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Abstract 



A measurement of the inclusive WW+WZ diboson production cross section in proton- 
proton collisions is reported, based on events containing a leptonically decaying W 
boson and exactly two jets. The data sample, collected at = 7TeV with the 
CMS detector at the LHC, corresponds to an integrated luminosity of 5.0 fb -1 . The 
measured value of the sum of the inclusive WW and WZ cross sections is c(pp — > 
WW + WZ) = 68.9 ± 8.7 (stat.) ± 9.7 (syst.) ±1.5 (lum.) pb, consistent with the stan- 
dard model prediction of 65.6 ± 2.2 pb. This is the first measurement of WW+WZ 
production in pp collisions using this signature. No evidence for anomalous triple 
gauge couplings is found and upper limits are set on their magnitudes. 
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The gauge symmetry of the standard model (SM) fixes the triple gauge boson couplings that 
determine the self-interactions of W and Z bosons. The pair production of vector gauge bosons 
allows a direct test of the electro weak sector of the SM [lj. Observation of anomalous triple 
gauge boson couplings would be an indication of physics beyond the SM. 

We report the first measurement of WW+ WZ diboson production in pp collisions in the semilep- 
tonic final state at the Large Hadron Collider (LHC), where one W boson decays leptonically 
(£v, with I = e, ji) while the other boson (W or Z) decays hadronically (jj), giving rise to two en- 
ergetic jets in the final state. Previous measurements in this channel at the Tevatron pp collider 
include the recent CDF [2] and DO ||3j HJ results. The advantage of reconstructing WW+WZ in 
the tvjj decay mode over the purely leptonic final states |5rEl is the larger branching fraction of 
W and Z bosons to quarks. This advantage is partially offset by the larger backgrounds in the 
ivjj channel, coming mainly from W+jets production. In contrast to the fully leptonic decay 
of WW pairs, the semileptonic process permits a direct measurement of the boson transverse 
momentum (px)- The sensitivity of WW production to the WW7 coupling and of WW and WZ 
production at high boson transverse momentum to the WWZ coupling makes these processes 
particularly useful as a probe of anomalous triple gauge boson couplings. 

The data correspond to an integrated luminosity of 5.0 ± 0.1 fb _1 collected in 2010 and 2011 
with the Compact Muon Solenoid (CMS) detector in pp collisions at y/s = 7 TeV at the CERN 
LHC. The CMS experiment |9] uses a right-handed coordinate system, with the origin at the 
nominal interaction point, the X axis pointing to the center of the LHC ring, the y axis pointing 
up, perpendicular to the plane of the LHC ring, and the z axis along the counterclockwise 
beam direction. The polar angle 9 is measured from the positive z axis and the azimuthal angle 
(p is measured in the x-y plane. The pseudorapidity is defined as rj = — ln[tan(0/2)]. The 
central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diameter, 
providing a magnetic field of 3.8 T. Within the field volume are silicon pixel and strip trackers 
and several calorimeters. The tracking system covers the range \t]\ < 2.5 and provides a track 
momentum resolution of 1% at 100 GeV. The lead tungstate crystal electromagnetic calorimeter 
(ECAL) covers \rj\ < 3 with an energy resolution of about 3%/ \/ E, where E is in GeV IITOl . The 
brass /scintillator hadron calorimeter (HCAL) covers \t]\ < 3.0 with an energy resolution of 
100%/ VE. The muon system consists of gas-ionization detectors inside and around the steel 
return yoke, and is capable of reconstructing and identifying muons within \ t]\ < 2.4. Extensive 
forward calorimetry complements the coverage provided by the barrel and endcap detectors. 
The CMS detector is nearly hermetic, allowing for measurements of the missing transverse 
energy (E™ ss ) in the event. A two-tier trigger system selects the events of interest. 

The data were collected with a suite of single-lepton triggers mostly using pj thresholds of 

24 GeV for muons and 25-32 GeV for electrons. To preferentially select events with on-shell W 
bosons, the single-electron triggers also require minimum thresholds on E™ ss in the range 0- 

25 GeV and on the transverse mass mj of the electron plus Eij? iss system in the range 0-50 GeV. 
The overall trigger efficiency is about 94% (90%) for muon (electron) data, with a small depen- 
dence (a few percent) on pj and rj. Simulated events are corrected for the trigger efficiency as 
a function of lepton pj and rj, and in the case of electrons also as a function of E™ 1SS . Simulated 
events are used to develop and validate the methods used in the analysis. 

The MadGraph5 1.3.30 [1 1 J event generator produces parton-level events with a Wboson and 
up to four partons on the basis of matrix-element (ME) calculations. The ME-parton shower 
(ME-PS) matching scale ]i is taken to be 20 GeV |T2"fl , and the factorization and renormalization 
scales are both set to q 2 = M 2 ^ + p TW . Samples of tt and Drell-Yan events are also gener- 
ated with MadGraph. Single-top production is modeled with POWHEG 1.0 (l3|. Multijet and 
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diboson samples (WW, WZ, ZZ) are generated with PYTHIA 6.422 JTU. PYTHIA provides the 
parton shower simulation in all cases, with parameters of the underlying event set to the Z2 
tune [15 j. The set of parton distribution functions used is CTEQ6LL Jl6|. A GEANT4-based 
simulation |T7| of the CMS detector is used in the production of all Monte Carlo (MC) sam- 
ples. Multiple proton-proton interactions within a bunch crossing (pileup) are simulated, and 
the triggers are emulated. All simulated events are reconstructed and analyzed as measured 
collision events. 

Events are selected with one well-identified and isolated lepton (muon or electron), large miss- 
ing transverse energy, and exactly two high-px jets. Muons are reconstructed within \t]\ < 2.1 
with the inner tracker and the muon system [fT8|. Electrons are reconstructed from tracks in 
the tracker pointing to energy clusters in the ECAL, within \tj\ < 2.5, excluding the transition 
region between the barrel and endcap, 1.44 < \tj\ < 1.57 p9j. Muons (electrons) are required to 
have px greater than 25 GeV (35 GeV). The lepton candidates are required to be consistent with 
having originated from the primary vertex of the event, which is chosen to be the vertex with 
the highest p | of its associated tracks. According to the simulation, this requirement provides 
the correct assignment for the primary vertex in more than 99% of the cases in both signal and 
background events. Charged leptons from W boson decays are expected to be isolated from 
other activity in the event. The sum of transverse momentum or energy in the tracker, ECAL, 
and HCAL, within a surrounding cone of AR = (Ar]) 2 + (A</>) 2 < 0.3, excluding the lepton 
itself, is required to be less than 10% of the measured pj of the muon, or less than 5% of the 
measured pj of the electron. Here Arj and A(p are the differences in pseudorapidity and in 
azimuthal angle, respectively. To reduce the backgrounds from fully leptonic decays, such as 
Drell-Yan and electroweak diboson processes, we exclude events in which there is any other 
loosely identified lepton (with pj > 10 GeV for muons and pj > 20 GeV for electrons) in the 
event. 

Jets are reconstructed from calorimeter and tracker information using a particle-flow tech- 
nique that combines information from several subdetectors |20|. The anti-fcx clustering algo- 
rithm plll22| with a distance parameter of 0.5 is used. Jets that overlap with isolated leptons 
within AR = 0.3 are not considered. Jet-energy corrections are applied to account for the non- 
linear energy response of the calorimeters and for other instrumental effects 11231 . These correc- 
tions are based on in situ measurements using dijet, 7+jet, and Z+jet data samples |24l| . Pileup 
collisions and the underlying event add to the energy of the reconstructed jets. The median en- 
ergy density from pileup is evaluated in each event and the corresponding energy is subtracted 
from each jet [25]. In addition, charged tracks that do not originate from the primary vertex are 
not considered for jet clustering II26II . We verified that these procedures sucessfully remove the 
dependence of jet response on the number of interactions in a single event. A jet-quality re- 
quirement, primarily based on the energy balance between charged and neutral hadrons in a 
jet, is applied to remove poorly reconstructed jets. Only events containing exactly two jets with 
pj > 35 GeV and within \t]\ < 2.4 are selected for the analysis. To reduce contamination from 
tt background, events are discarded if one or more jets pass high-efficiency b-quark jet iden- 
tification criteria based on the presence of a secondary vertex within the jet [27j. An accurate 
E™ 1SS measurement is essential to distinguish the W signal from multijet backgrounds and to 
reconstruct the full event kinematics of the WW system. We use E™ 1SS measured in the event 
using the full particle-flow reconstruction |28f and require E™ ss > 25 (30) GeV for the muon 
(electron) channel. To reduce the background from processes that do not contain W — > iv de- 
cays, we require that the transverse mass of the W candidate exceed 30 GeV (50 GeV) in muon 
(electron) data 1291 . 

We measure the dijet mass {tnu) distribution, as shown in Figure [TTa). The relative contribu- 
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Table 1: Treatment of background nijj shapes and normalizations in a fit to the data. The cross 
section values are calculated with the programs cited on the corresponding rows. The back- 
ground normalizations are constrained to Gaussian distributions with the listed central values 
and widths. The treatment of multijet events is described in the text. 



Process 


Shape 


Constraint on normalization 


Diboson (WW+WZ) 


MC 


Unconstrained 


W+jets 


MC 


31.3nb±5%(NLO) El 


tt 


MC 


163pb±7% (NLO) |3ll 


Single top 


MC 


85 pb ± 5%(NNLL) BHU 


Drell-Yan+jets 


MC 


3.05 nb ± 4.3% (NNLO) flU 


Multijet (QCD) 


data 


E™ iss fit in data 



tions of the known SM processes are determined using an unbinned maximum-likelihood fit 
over the mass range 40-150 GeV. The fit is performed separately for the muon and electron 
channels since their background compositions differ. Table [TJlists the SM processes included in 
the fit. The normalization of the diboson WW+WZ contribution is a free parameter. The nor- 
malizations of the background components are allowed to vary within Gaussian constraints 
around their central values. For multijet events, this central value is obtained from an indepen- 
dent two-component fit to the E™ 1SS distribution which determines the corresponding fraction 
in the data Il29l . The fit uncertainty is used as a constraint on the multijet contribution. The 
central values for all other processes are obtained from next-to-leading-order (NLO) or higher- 
order calculations, and the constraints are taken from the theoretical uncertainties. With the 
exception of multijet production, the shape of the nijj distribution for all processes is obtained 
from simulation. Multijet events contribute to the total background when jets are misidentified 
as isolated leptons. Their nijj shape can be derived from data events with lepton candidates 
that fail the isolation requirements. 

The nijj spectrum of the dominant W+jets component is described using the shape from MAD- 
Graph simulation after taking into account the uncertainties due to the factorization and 
renormalization scale (both equal to q) and ME-PS matching scale \i |36|: 

fW+jets = * -^W+jets (Fo^' 2 ) + $ -^W+jets (f*' 2 ' *7o) 

+ (l-a-/J) JW+jets(> 2 /<7o), (!) 

where ^w+jets denotes the m» shape from simulation. The parameters y.Q and qo (q') cor- 
respond to the default (alternative) values of ji and q, respectively. The parameters a and 
jS are free to vary during the fit and remain within the physical ranges (0 < a, f> < 1 and 
1 — ol — /3 > 0). We take ]i' = 2}Iq or 0.5^o (<?' = 2^o or 0.5^o)/ depending on which alternative 
sample provides a better fit to the data. Thus, the fit probes variations of a factor of two in both 
f/ and q. 

Figure [T|a) shows the observed nijj distribution for both channels combined, together with the 
fitted projections of the contribution of various SM processes. FigureJTJb) shows the same distri- 
bution after subtracting all SM contributions from data except for WW+WZ events. Figure [TJc) 
shows the pull distribution, i.e. the normalized residual defined as (data — fit) / (fit uncertainty), 
where the fit uncertainty is computed at each data point by propagating the uncertainty in the 
normalization coefficients. The yields of various SM components, as determined by the fit, are 
reported in Table [2] 

We validated the fit procedure by performing pseudo-experiments. In each experiment, we 
generate the nijj pseudo-data for the SM processes, taking into account the correlations between 
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Figure 1: (a) Distribution of the dijet invariant mass in data, with the binning chosen based on 
the resolution and fit projections of the relevant components overlaid, (b) The dijet invariant 
mass after subtraction of all components except the electroweak WW+WZ processes. The er- 
ror bars represent the statistical uncertainties and the hatched bands represent the systematic 
uncertainties, (c) The normalized residual or pull: (data — fit)/ (fit uncertainty). 

the yields, and then perform a fit to each pseudo-data m« distribution. The results for both the 
muon and electron channels indicate that there is a small bias (-8.6% and -6.6%) in the WW+WZ 
yield, corresponding to less than 0.4 standard deviations, and that the fit slightly overestimates 
the uncertainty on the yield. These effects are corrected for in the final result. The validation 
procedure shows that biases in all background yields and errors are small. The fit results for the 
background components are statistically consistent with the expectations, with the exception of 
W+jets, where 12% fewer events than expected are observed. Overall, the approach produces 
a high quality model of the data (Fig. [l|a)), where the pull distribution is consistent with 
(Fig.JTJc)), and allows us to extract the diboson peak (Fig.JTpb)). 

Systematic uncertainties arising from the jet energy are estimated from W bosons decaying 
hadronically in a sample of semileptonic tf events. The mean and resolution of the recon- 
structed dijet mass distribution in data agree to within 0.6% of the expectations from simula- 
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Table 2: Event yields determined from a maximum-likelihood fit to the data. The total uncer- 
tainty is computed using the full covariance matrix. Owing to a higher kinematic threshold the 
product of acceptance x efficiency is smaller for the electron channel. The term Ae includes W 
and Z branching fractions J37|. 



Process 


Muon channel 


Electron channel 


Diboson (WW+WZ) 


1900 ±370 


800 ± 310 


W plus jets 


67380 ± 590 


31640 ± 850 


tf 


1660 ± 120 


950 ± 70 


Single top 


650 ± 30 


310 ± 20 


Drell-Yan+jets 


3610 ±160 


1410 ±60 


Multijet (QCD) 


300 ± 320 


4190 ±870 


Data 


75419 


39365 


Fitx 2 /N dof (probability) 


9.73/12 (0.64) 


5.30/12 (0.95) 


Acceptance x efficiency (Ae) 


(5.15 ±0.24) x 10~ 3 


(2.63 ±0.12) x 10~ 3 


Expected WW+WZ yield from simulation 


1700 ± 60 


870 ± 30 



tion. A small difference in E™ 1SS resolution II28II between data and simulation affects the signal 
acceptance at the 0.5% level. Further systematic uncertainties on the signal yield are due to the 
uncertainty on the trigger efficiency in data (1%), and on the lepton reconstruction and selection 
efficiencies (2%) [29 1. The uncertainty due to the b-jet veto is negligible. The uncertainty in the 
luminosity measurement is 2.2% [38]. The uncertainty in acceptance arising from theoretical 
uncertainties, including parton distribution functions and additional jet rejection, is 4%. 

As listed in Table [2} we observe 2700 ± 340 (stat.) ± 360 (syst.) WW+WZ events, in agreement 
with the SM expectation. This result corresponds to a significance of 8.8 standard deviations 
when computed using a simple likelihood ratio, where the nuisance parameters [39J are fixed 
to their fitted values. Using the profile likelihood ratio |39| , where the nuisance parameters are 
allowed to vary, the significance becomes 4.3 standard deviations. We compute the WW+WZ 
total cross section as a = Nsi g / (AeC), where Nsig is the number of extracted signal events, A 
is the signal acceptance corrected for the branching fractions, e is the overall efficiency for event 
selection, and C is the integrated luminosity. In the acceptance calculation we assume the SM 
value for the WW to WZ production ratio. The values of N$i g and Ae are given in Table [^sepa- 
rately for the muon and electron channels. Combining the results from the muon and electron 
channels, we obtain a(pp -> WW ± WZ) = 68.9 ± 8.7 (stat.) ± 9.7 (syst.) ± 1.5 (lum.) pb, which 
is in agreement with the NLO prediction [40J, 65.6 ± 2.2 pb, that includes the contribution from 
gg — > WW. The total cross section computed in the muon channel, 73.8 ± 15.1 pb, is consistent 
with that obtained in the electron channel, 60.8 ± 21.5 pb, where the statistical and systematic 
uncertainties have been combined in quadrature. 

Measurements of electroweak diboson production can be translated into measurements of 
gauge boson self-couplings, which are among the most fundamental aspects of the SM. At 
the leading order, only s-channel qq annihilation diagrams have a three-boson vertex involv- 
ing WW7 and WWZ couplings in WW production, and WWZ coupling in WZ production. 
Physics beyond the SM can modify these couplings, leading to observable differences in the 
cross section and the kinematic distributions fl4"T| . We search for anomalous triple gauge cou- 
plings (aTGCs) using an effective Lagrangian described by the following HISZ (Hagiwara, 
Ishihara, Szalapski, Zeppenfeld) parametrization without form factors [lj: A 7 = Az = A, 
Akz = Agf — Ak 7 • tan 2 0w. We use the dijet pj distribution, shown in Fig. [2J as the observable 
after requiring 75 GeV < m,-; < 95 GeV to enhance signal purity. The dependence of the dis- 
tribution on specific aTGCs is modeled by reweighting the pythia simulation of WW+WZ to 
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Figure 2: Dijet pj distributions for (a) the muon and (b) the electron channels after full selec- 
tion and with the requirement 75 GeV < rrijj < 95 GeV. The stacked histogram shapes are 
taken from simulation or, where applicable, from data-driven estimates. They are normalized 
according to the fit to the observed ntu spectrum in data. Below we show the Data/MC ratio 
with the (dashed) red lines corresponding to the shape uncertainty. The last bin includes the 
overflow. 



the MCFM [30 1 NLO predictions. We account for systematic uncertainties arising from luminos- 
ity, signal selection efficiency, signal shape, and from the normalization and shape of the SM 
processes. We find no evidence for aTGCs. Given the tight bound on the parameter Agf |37|, 
we assume the SM value (Agf = 0) and set limits on the two parameters A and Ak t Exclu- 
sion limits at 95% confidence level (CL) in the two-dimensional A-Ak 7 plane, computed using 
the modified frequentist CL$ l|39l l42l technique with profile likelihood as the test statistic, are 
shown in Fig. [3] We obtain the following one-dimensional observed 95% CL limits assuming 
the SM value for the other parameter: -0.038 < A < 0.030, -0.11 < Ak 7 < 0.14. These limits 
are competitive with, and in some cases improve upon, the sensitivity of the combined LEP ex- 
periments listed in Refs. 11371 |43T(46| . The ATLAS Collaboration recently reported limits in the 
fully leptonic channel for WZ [7] and WW [8J production. Limits obtained from fully leptonic 
channels are weaker due to the smaller branching fractions. 

In summary, a measurement of the sum of the inclusive WW and WZ production cross sections 
has been performed using events containing a leptonically decaying W and two jets. The mea- 
sured value for the cross section is 0"(pp — >• WW + WZ) = 68.9 ± 8.7(stat.) ± 9.7(syst.) ± 
1.5 (lum.) pb, which is consistent with the SM prediction. This is the first measurement of 
WW+WZ production in pp collisions using this signature. No evidence for anomalous triple 
gauge couplings is found, and limits are set on their magnitudes. 

We congratulate our colleagues in the CERN accelerator departments for the excellent perfor- 
mance of the LHC and thank the technical and administrative staffs at CERN and at other CMS 
institutes for their contributions to the success of the CMS effort. In addition, we gratefully 
acknowledge the computing centers and personnel of the Worldwide LHC Computing Grid 
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by the solid dot. 
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